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We discuss various aspects of Conformal Field Theories on the Lattice. We mainly in- 
vestigate the SU (3) gauge theory with Nf degenerate fermions in the fundamental rep- 
resentation, employing the one-plaquette gauge action and the Wilson fermion action. 

First we make a brief review of our previous works on the phase structure of lattice 
gauge theories in terms of the gauge coupling constant and the quark mass. We thereby 
clarify the reason why we conjecture that the conformal window is 7 < Nf < 16. 

Secondly, we introduce a new concept, "conformal theories with IR cutoff' and point 
out that any numerical simulation on a lattice is bounded by an IR cutoff Aju. Then 
we make predictions that when Nf is within the conformal window, the propagator of a 
meson G{t) behaves at large t, as G(t) = c exp {—m[jt)/t°', that is, a modified Yukawa- 
type damping form, instead of the usual exponential damping form exp (— mfft), in the 
small quark mass region where mjj < cAjn is satisfied: mj{ is a generic mass of a hadron 
such as the pion mass and c is a constant of order 1. This holds on an any lattice for any 
coupling constant g, as far as g is between and g*, where g* is the IR fixed point. We 
verify that numerical results really satisfy the predictions for the Nf = 7 case and the 
Nf = 16 case. 

Thirdly, we discuss small number of flavors (Nf = 2 ~ 6) QCD at finite tempera- 
tures. We point out theoretically and verify numerically that the correlation functions 
at T/Tc > 1 exhibit the characteristics of the conformal function with IR, cutoff, an 
exponential damping with power correction. 

Investigating our numerical data by a new method which we call the "micro-analysis" 
of propagators, we observe that the Nf = 7 case and the Nf = 2 at T ~ 2 Tc case are 
similar to each other, while the Nf = 16 case and the Nf = 2 at T = 10^ ~ lO^Tc cases 
are similar to each other. 

Further, we observe our data are consistent with the picture that the Nf = 7 case 
and the Nf = 2 at T ~ 2Tc case are close to the meson unparticle model. On the 
other hand, the Nf = 16 case and the A?^ = 2 at T = 10^ ~ lO'^ Tc cases are close 
to the fermion unparticle model. All results are consistent with naive physical intuition 
and solve long standing issues at high temperatures such as why the free energy at high 
temperatures does not reach the Stefan-Boltzmann ideal gas limit even at T = 100 Tc. 

Keywords: Conformal Theories, Conformal Window, IR cutoff, Unparticle 
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1. Introduction and Summary 

Recently there has been a lot of interest in finding a conformal or near-conformal the- 
ory, which could play a key role in constructing a theory Beyond the Standard Model 
(BSM) of particle physics.^ Non-perturrbative investigation of non-supersymmetric 
conformal theories in four dimensional space-time is also ardently desired. 

We already described in the abstract partly the objectives, the target, the 
method and the results of this work. In order to avoid to repeat the same things, we 
will skip what was described in the abstract and will describe remaining important 
matters. 

The strategy we take is as follows: We define continuum conformal theories as 
the continuum limit of lattice theories. This is a constructive field theory approach. 
First, we propose the properties of continuum conformal theories from knowledge 
of the continuum super-symmetric conformal theories and RG argument .1^ Then, 
we conjecture the properties of lattice theories. 

Lattice theories are defined on a Euclidian four-dimensional lattice with lattice 
size Nx, Ny, Nz, Nt and lattice spacing a. For simplicity, we assume = Ny = 
Nz = N and A''^ = rN, with r a constant. The continuum limit of the theory is 
defined by a — > with iV — ^ oo keeping L — N a constant. When L = finite, the 
continuum limit defines a theory with an IR cutoff. 

We propose that the propagator of a meson G{t) behaves at large t, as G{t) = 
c exp (— r7iijt)/t" for a small enough quark mass. We note that this form of the 
propagator in the coordinate representation corresponds to a cut instead of a pole 
in the momentum representation. 

After we have verified that numerical results satisfy the predictions for the Nf = 
7 case and the Nf ~ 16 case, we propose a new method of analysis of propagators 
which we call the " micro-analysis of propagator" : Writing the propagator of a meson 
as G{t) = c exp {—m{t)t)/t"^^\ investigate the t dependence of m{t) and a{t). The 
m{t) and a{t) are reflected by the evolution of RG transformation flow and therefore 
contain useful information of the dynamics. 

Carefully investigating the m(t) and a(t), we observe that our data are consistent 
with the picture that the Nf = 7 theory is close to the meson unparticle model, 
while the Nf = 16 theory is close to the fermion unparticle model. We also observe 
the similarity between the large Nf case and the small Nf at T > Tc case. 

2. Stage and Tools 

We employ the Wilson quark action and the standard one-plaquette gauge action. 
For fermion field an anti-periodic boundary condition is imposed in the time direc- 
tion, otherwise periodic boundary conditions. The theory is defined by two param- 
eters; the bare coupling constant go and the bare quark mass mo. 

We usually use, instead of the bare gauge coupling constant go and the bare 
quark mass toq, P defined by /3 = 6/gg and K — 1/2 (moa -I- 4), which is called the 
hopping parameter. 
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In order to investigate properties of the theory in the continuum limit, one has 
to first clarify the phase structure of lattice QCD at zero temperature to identify 
a UV fixed point and/or an IR fixed point. When Np < 16, the point go = and 
mo = is a UV fixed point. Therefore a theory governed by this fixed point is an 
asymptotically free theory. We restrict ourselves to the theory defined by this UV 
fixed point in this articles. 

When there is an IR fixed point^J at finite coupling constant g on the massless 
line which starts from the UV fixed point, the long distance behavior is determined 
by the IR fixed point. This defines a conformal theory. The minimum number of 
flavors Nf for the existence of such an IR fixed point is denoted by N^. Usually 

< ^ 16 is called the conformal window. 

The local meson operator is denoted by Gnit)- When the theory is in the con- 
fining phase, it decays exponentially at large t 

Gnit) = c exp{-mHt) (1) 

for a free boundary, where mn is the mass of the ground state hadron contained in 
the channel H. 

The boundary conditions we take imply 

G{t) = c (exp(-TOH<) + exp(-mff (L( - t))) (2) 

for t — Lt/2. From this behavior we extract the energy of the ground state mn in 
the channel H. 

The quark mass is defined through Ward-Takahashi identities 



_ (OIV4A4IPS) 
^" 2(0|P|PS) ' ^ ' 

where P is the pseudo-scalar density and ^4 the fourth component of the local axial 
vector current. Renormalization constants are suppressed, since this expression is 
enough for our purpose. 



3. Phase Diagram: a brief survey of our previous v^rorks 

In our earlier article,'^ we investigated the phase structure of the SU(3) and SU(2) 
gauge theory with various numbers of flavors Nf of Wilson fermions in the funda- 
mental representation, and conjectured that the conformal window is 7 < Np < 16 
for the SU(3) and 3 < Np < 10 for the SU(2), based on numerical results for various 
numbers of flavors. 

Although we would like to extract the phase diagram at infinite volume and at 
zero temperature, we have to make simulations at finite volume N^ and at finite 
Nf — rN . Therefore the phase diagram is a three dimensional space parameterized 
by, go, ruQ and N. Finally we have to take the limit N infinity. 
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3.1. Nf <6 

When TV is finite, the phase diagram is given as in Fig.IIIa). In 1996, we studiecP 
the finite temperature transition of QCD for Np —2, 3, 6. Our results implies that 
the transition in the chiral limit is of first order for Nf — 3 and 6. The order of the 
chiral transition in the Nf = 2 case is more subtle. We will discuss this issue later. 

At K = which corresponds to the pure gauge theory, the phase transition is of 
first order. Thereby the transition line connects between the critical point at K = 
and the chiral transition point at K = Kc- 

As lattice size N increases, the phase transition points both at quenched QCD 
and at massless QCD move towards larger (3 as shown in Figlljb). Eventually in 
the limit — )■ oo, all the region is covered with the confinement phase as shown in 
FigEIa). 



3.2. 7 <Nf <1Q 

When N is finite, we found that the phase diagram is given as in Fig. 2(a). Due to 
the lack of chiral symmetry the structure of phase diagram is complicate. 
The salient facts are the following: 

First, the massless line which starts the UV fixed point moves towards smaller 
/3 and hits the bulk transition. The massless line belongs to a deconfining phase all 
through the region. 

Secondly, in the confinement phase which exists at small /3 region there is no 
massless line. When we decrease the quark mass from the quenched QCD, we en- 
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Fig. 1. Phase diagram for Nf < 6. (a) For not so largeTV. (b) For larger A'^. 
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Fig. 2. Phase diagram for 7 < Nf < 16. (a) For not so largeAf. (b) For larger A'^. 
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counter a bulk transition at finite quark mass. Thus this confinement phase is irrel- 
evant to the continuum theory. 

As lattice size N increases the phase diagram changes as shown in Figl^l^b) and 
FiglHl^b). It should be noted that there is still no massless line which belongs to 
a confinement phase. This implies that even in the continuum limit, the massless 
line does not belong to the confinement phase. On the other hand, if there would 
be no IR fixed point, the massless line should belong to a confining phase in the 
continuum limit. From this reasoning we conclude that the Nf — 7 theory is a 
conformal theory. Thus we conjecture that the conformal window is 7 < Nf < 16. 

4. Strategy for Identification of tiie Conformal Window 

The logics in the previous work is an indirect way to identify the conformal window. 
Certainly more direct ways are desirable. The most direct and solid way is to identify 
the IR fixed point. However, for small iV^ the IR fixed point exists deep in strong 
coupling region, and it becomes difficult to identify the IR fixed point. 

An alternative way is to find directly the characteristics of conformal theories. 
This is the way we take in this work. 

5. Conformal Field Theories 

In this section we address ourselves to a basic and fundamental issue of what hap- 
pens when the number of flavor Nf is within the conformal window. 

The strategy we take is as follows: We define continuum conformal theories as 
the continuum limit of lattice theories. First, we propose the structure of continuum 
conformal theories from knowledge of the continuum super-symmetric conformal 
theories and RG argument. Then, we conjecture the structure of conformal theories 
on the lattice. 

5.1. Continuum Conformal Field Theories 

We discuss separately the L — oo case and the L = finite case, since there exists a 
clear difference between the two cases. 
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Fig. 3. Phase diagra for large N. (a) For Nj < 6. (b) For 7 < Nf < 16. 
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5.1.1. The case L — oo 

This defines a continuum theory on i?** Euclidean four-dimensional space, which 
corresponds to the flat infinite space at zero temperature. In this case the IR cutoff 
A/fl is zero. 

Therefore the relevant quantity is only the running coupling constant and 
there are no physical quantities which have physical dimensions. 

The "conformal region" exists only on the massless line. All other region is the 
"confinement region". The phase structure is shown in Fig|4lja). In super-symmetric 
theories this kind phase structure is known. The shaded region in the strong coupling 
region does not exist in the /3 — ruq plane. When the phase structure is described in 
terms of 13— K, the corresponding phase belongs to a region for Wilson doubles (Aoki 
phase). Therefore when it is mapped to the /? — planes, the region corresponding 
to the shaded one does not exist. 

The propagator on the massless quark line takes the form 



(4) 



since, as mentioned, there is no physical quantities with physical dimensions. 
If we take the coupling constant g = g* at UV scale ji, 

a(t) = 3-27*, (5) 

with 7* being the anomalous mass dimension ^ at g = g* . The theory is scale 
invariant (and shown to be conformal invariant within perturbation theory.l^ See 
also e.g. Nakayama^ and references therein from AdS/CFT approach). In general, 
the power a varies with t: 

a = a{t). (6) 

When t is small, the correlation function is governed by the UV fixed point. 
There the "meson" state is a free quark and an anti-quark state. Thus the power is 
given by 

a(t) = 3. (7) 
When t is large, it is governed by the IR fixed point. The exponent is given by 

a{t) = 3 - 27*. (8) 



IR 
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confinement 




P = 00 |3 = o p : 

(a) (b) 

Fig. 4. Phase diagra for large = 00. (a) Aju = 0. (b) Ajjj = finite 
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The exponent of the power a in IR is given by the renormaUzation equation or 
derived just by dimension counting. 

5.1.2. The case L — finite 

This defines a continuum theory on T'^ Euchdean four-dimensional space, which 
corresponds to a compact space at finite temperature. In this case the IR cutoff 
Ajji = 1/L is finite. 

The renormaUzation transformation to IR does stop to evolve at the scale of IR 
cutoff A/jj. Thus the relevant physical quantities we may choose are 

AcFT,-AlR,mH (9) 

Here Ac ft is a scale parameter for the transition region from UV to IR, and tuh 
is a generic mass of a meson such as the pion and the p meson. 

Now we make a proposal. The propagator of a meson G{t) behaves at large t, as 

Git)^c^^^, (10) 
that is, a modified Yukawa-type damping form in the region 

rriH < cAiR, (11) 
instead of the usual behavior of the propagator in the " confinement region" , 

G{t) = ccxp{-mt). (12) 

Here A/^ is an IR cutoff, rriH is a generic mass of a hadron such as the pion 
mass., and c is a constant of order 1. The "conformal region" is depicted in FigHJb). 
The relation ea.([TT|) reduces to 

<c(l/L)(i+^*), (13) 

using the relatiorP 

mjj = cmj/(i+'^*). (14) 

The phase structure with L = finite smoothly approaches the case L = oo, as L 
varies to oo. The form of the propagator also smoothly becomes 

G(t) = cl, (15) 

in the massless limit. 

As is the case in the Am — case, the power a varies with t: 

a = a{t). (16) 

When t is small, 

a{t) = 3. (17) 
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On the other hand, when t is large {t ^ 1/Ac_ft), the behavior of a{t) is 
comphcate: When m is small enough and if mi ^ 1, the power is given by 

a(t) = 3-27*. (18) 

If we take larger t such that mt ^ 1, the a{t) takes a value which depends on the 
dynamics of the theory. We give explicit examples when we discuss some models. 

The concept of " conformal theory with IR cutoff" and " conformal region" de- 
termined by eqations (fTO|) and (jlip are our key proposals. 

5.2. Conformal Field Theories on the Lattice 

We are able to perform numerical simulations only on a finite size lattice. Thus the 
system is automatically bounded by the IR cutoff A//j ^ l/{Na). In this case, the 
relevant physical parameters we may use are Ac ft, A/^ and ttih- 

Our main target is first to verify the transition of a meson propagator from 
an exponential damping form to a modified Yukawa-type, that is, an exponential 
form with power correction given by eq. (|10|) . at the critical hadron mass given by 
eq. ini). 

When the lattice size A^i increased to iV2 = s iVi with /3 being fixed, then 

mg;'*'™'(/3,7V2) = l/sTO^***^'''(/3,iVi). (19) 

The critical hadron mass decreases. Thus we have to carefully choose the simu- 
lation parameters in order to find the "conformal region" and the Yukawa- type 
propagators. 

6. Numerical Simulations 

6.1. Algorithm and Parameters 



Table 1. Parameters of Simulations 



Nf /3 
7 6.0 



K 


mq 


Nf 13 
16 11.5 


K 


mq 


Nf 
2 


6.5 


K 


mq 


0.1400 


0.22 




0.126 


0.22 






0.146 


0.035 


0.1446 


0.084 




0.127 


0.19 




10.0 






0.1452 


0.062 




0.130 


0.10 






0.135 


0.028 


0.1459 


0.045 




0.1315 


0.055 




15.0 






0.1472 


0.006 




0.13322 


0.003 






0.130 


0.046 



We have made simulations mainly for the cases Nf — 7 and Nf — 16, which are, 
we think, the boundaries of the conformal window. The algorithm we employ is the 
blocked HMC algorithnP for Nf ^ 2N and the RHMC algorithnpl for Nf = 1, 
in the case Nf = 2 N + 1. The parameters we have simulated are listed in Tablel. 
(The parameters for Nf = 2 at high temperatures are also listed.) We have chosen 
the hopping parameter K in such a way that the quark mass takes the value from 
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0.25 to 0.0. It might be worthwhile to mention that we are able to take K which 
corresponds to a negative quark mass crossing the zero quark mass, since the phase 
is a deconfining phase and therefore the pion mass is not zero there. The coupling 
constant of the IR fixed point for Nf = 16 is /3 = 11.48 in one-loop approximation. 
This is the reason why we choose fi = 11.5 for Nf = 16. 

The lattice sizes we choose are 8^ x 32, 16^^ x 64 and 24 x 96. The run-parameters 
are chosen such that the acceptance of the global metropolis test is about 70%. 
The statistics is in general 1,000 MD trajectories for thermalization and 1,000 MD 
trajectories for measurement. 

7. Numerical Results 

Let us first discuss the results for Nf = 7 case on a 16"^ x 64 lattice. FiglS] shows the 
effective mass plots for pions with the local-sink local-source (black squares), local- 
sink doubly-smeared-source (red circles) and local-sink doubly-wall-source (green 
triangles) cases. The effective masses in the relatively large quark mass case shown 
on the left panel, ruq = Q.2^{K = 0.1400), clearly show the plateau at t = 22 31. 
On the other hand, in the = 0.045 {K = 0.1459) case shown on the right panel, 
there is no plateau up to i = 31. The effective mass is slowly decreasing. 

The effective masses for all cases with iriq < 0.084 {K > 0.1446) exhibit similar 
behavior to K = 0.1459. We have conjectured that propagators decay exponentially 
with power correction. Thus, to verify the conjecture, we fit the local-local data 
with the Yukawa-type for the fitting range t = [15 : 31]. The fits for K — 0.1459 
and K = 0.1452 are depicted in Figl6l The both fits well reproduce the data with 
range t = [15 : 31] with a = 0.8(1). 

In the case of Nf = 16, also, the effective masses on a 16'^ x 64 lattice exhibit 
the modified Yukawa-type for all cases when K > 0.130. The results of the fit to the 
local-local data for K = 0.130 and K = 0.1315 with the Yukawa-type for the fitting 
range t = [15 : 31] are plotted in FiglTl The fits well reproduce the data. However, 
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Fig. 5. The asymptotic behavior of propagators in the Nf = 7 case on a 16^ X 64 lattice at 
/3 = 6.0. (a) Exponential damping form.(b) Modified Yukawa-tye damping form. 
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if we carefully compare the fits of Nf — 7 and Nf — 16 cases, the fit of Nf — 16 
case at large t is not so good as Nf = 7 case. We will discuss this point later. 

The data on a 24^ x 96 lattice for both cases of iV/ = 7 with K > 0.1459 and 
Nf = 16 with K > 0.130 exhibit the modified Yukawa type decay. 

Thus we have verified that the propagators decay with the Yukawa-type damping 
form for the Nf — 7 and Nf = 16 cases, for the small quark mass. In term of the 
pion mass the critical mass can be expressed as mn < cAjji. Here mn is the pion 
mass and we take Air = {N^ x NtY/^. There is no a priori definition of Ajr in the 
case of asymmetric lattices. This is one of choices. With this definition, we estimate 
the constant c is c 11.0 for Nf — 7, and c ^ 12.4 for Nf — 16. 



K=0.1459, mq=0.045 K=0.1452, mq=0.062 

0.5 t V , 0.52 t 

■■, , loc(t)-loc(0) : loc(t)-loc(0) 
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0.45 - : 0.47 - %T,. 
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0.43 I - U 1 1 J 0,45 _ 

0.42 f ^ 0.44 f ' 

2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 

(a) (b) 

Fig. 6. Two examples of modified Yukawa-type damping form in the A'^^ = 7 case on a 16^ X 64 
lattice at 13 = 6.0 (a) K=0.1459. (b) K=0.1452. 
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Fig. 7. Two examples of modified Yukawa-type damping form in the Nf = 16 case on a 16^ X 64 
lattice at 13 = 11.5. (a) K=0.130. (b) K=0.1315. 
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8. Micro-analysis of Propagator 

As we have verified our conjecture tliat tlie propagators decay witli the Yukawa- 
type damping form for the small quark mass, the next issue we would like to solve 
is what kind of theory is defined. To do so, we propose a new method of analy- 
sis of propagators which we call the "micro- analysis" of propagators: Writing the 
propagator of a meson as 

G{t) = c exp {-m{t)t)/t"'^*^ (20) 

investigate the t dependence of m{t) and a{t). The m{t) and a{t) are reflected by 
the evolution of RG transformation flow and therefore contain useful information 
of the dynamics. 

We depict m{t) and a{t) in Figs. 9 ~ 19 for various cases. 

First, at small t (we disregard the a{t) and m{t) at t = 1 and 2, since they are 
affected by the boundary effect), 

a{t) ~ 3.0 (21) 

and 

m(t) ~ 2 mg (22) 
are well satisfied, when niq > 1/N. 

As t increases, m{t) increases from 2mq to mn- The behavior of a{t) is more 
complex and it reflects the dynamics. 

In the Nf = 7 case on the 16"^ x 64 lattice, the a{t) at large t is relatively stable 
and it seems that a{t) shows a plateau at t = 15 ~ 31, although it is noisy. The 
statics of the data on the 24 x 96 lattice is not enough to see the a{t), as it requires 
more statistics than the Yukawa-type fit with a wider range of t. We are storing 
configurations. 

The t dependence of a{t) in the Nf = 16 case is quite different from that of 
the Nf = 7 case. In the Nf = 16 case, although there is no plateau structure in 

a{t), there is a shoulder with a value around 1.5 at i = 12 ^ 24. on the 16^ x 64 
lattice. On the 24 x 96 lattice (the statistics is enough, since the Nf = 16 case is 
much faster than the Nf = 7 case), the t dependence of a{t) is different each other. 
In the case K = 0.130, the shoulder becomes very wide down to t ~ 32. In the 
case K — 0.1315 the shoulder is clearly seen for t = 12 ^ 25. on the other hand, 
the shoulder disappears and the plateau at a ~ 0.5 is seen for t > 24 in the case 
K = 0.13322. 

In order to find what kind of theory is defined, wc will investigate some models 
which we consider are good candidates to compare with the theory. To do so, we 
find it is better to discuss conformal theories at high temperatures before that. 
Therefore we examine conformal theories at high temperatures in the next section. 
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9. Conformal Field Theories at Finite Temperature 

Now let us consider lattice QCD at high temperatures for small Nf (2 < Nf < 6) 
where the chiral phase transition occurs at Tc- Note that the system is bounded 
by the infrared cutoff Am = l/T. The continuum limit of the theory at finite 
temperature is defined by a — > with N ^ oo and Nt —?' oo keeping Na = L 
and Nta = l/T constant. Taking the limit i — oo, the gauge theory at finite 
temperature T is defined. 

We first consider the case that the renormalized quark mass is zero. Then the 
renormalized coupling constant is the only relevant variable in the theory. A running 
coupling constant g(fi; T) at temperature T can be defined as in the case of T = 0. 
The following discussion can be applied to any method of the definition of the 
running coupling constant g{pL\ T). One way is the Wilson RG method. 

Although in several literatures (see e.g. ref. llip the running coupling constant 
g{lJ,;T) was discussed, it was mainly the UV behavior. Up to now it seems that it 
has not been discussed systematically the IR behavior. 

In the UV region the running coupling constant at finite T is essentially as the 
same as that of T = 0. Since the theory is asymptotically free, at small g region the 
following relation 



with c constant, is satisfied. The constant c depends on both of the schemes of 
g(/i;r = 0) and g(/i;T) in addition to Nf. It might be worthwhile to stress that 
the leading term is universal in the limit g 0. 

However, in the IR region g{^; T) may be quite different from ^(/i; T = 0), since 
the IR cut-off Km in the time direction is l/T. This is quite different from the 
T = case. 

Furthermore, the fact that when T/T^ > 1, the quark is not confined, implies 
that the running coupling constant g[^\T) cannot be arbitrarily large. This means 
that there is a fixed point in the fi function when T/Tc > I. 

Let us increase the temperature gradually from T = 0. 

As far as T < Tc, the (3 function is negative all through g. As the temperature is 
increased further, the form of j3 function will changes as in FiglHl (a) When T > Tc 



g^ifi; T) = g\ti- T = 0) + cg\^l■ T = 0), 



(23) 




(a) (b) 
Fig. 8. 13 function with T = constant, (a) T > Tc. (b) T > Tc. 
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but T ~ Tc, the jS function changes the sign from negative to positive at large g. 
(b) When T ^ Tc it changes the sign at smah g. 

The T/Tc dependence of the /? functions resembles the Nj dependence for large 
Nf. Nf = IQ^ T/Tc > 1 and Nj = 7 ^ T/Tc > 1. 

9.1. Numerical Results at High Temperatures 

We have made simulation in the Nf = 2 case on a 16'^ x 64 lattice for several sets 
of parameters. Let us discuss the three cases listed in Table. 1: /3 = 6.5 corresponds 
to T - 2Tc, /3 = 10.0 to T ~ 100 T^, and /S ^ 15.0 to T - lO^T^. 

We clearly see that the propagators in the three cases decay with the modi- 
fied Yukawa type damping form. Therefore we make the "micro- analysis" of the 
propagators. The m(t) and a{t) are depicted in Figs. [T4l [TSl and [TBI 

10. Correspondence between Nf and T 

We observe the t dependence of Tn{t) and a{t) in the case T ~ 2Tc are very similar 
to those in the Nf = 7 case. On the other hand, the t dependence in the T ^ lOOTc 
and T ~ lO^Tc cases are very similar to those in the Nf = 16 case. We plot in 
Figs. [171 HH and [191 the large Nf data on the right panel and the Nf = 2 at high 
temperature data on the left panel in order to compare them directly. 

The similarities between them are excellent. We could anticipate this from the 
similarities of the f3 functions: Nf = 16 ^ T/Tc > 1 and iV/ = 7 4^ T/T^ > 1. 

We point out a possibility that there is a correspondence between the large Nf 
with Aju and small Nf at high temperature T. Both of the large Nf with Am and 
the small Nf at high temperature T are described by two parameters in addition to 
N and Nt; P and vriq. In order to have a correspondence we have to map from one 
set of the parameters to the other set. Although It is too premature to discuss the 
details, the correspondence between Nf and T would be very intriguing and very 
useful. Although the flavor number Nf is a discrete integer, the temperature T is a 
continuous number. Thus the properties of large Nf could be deduced from those 
of small Nf at high temperatures with this correspondence. 

11. Models 

We investigate a few models to compare the results with the models. 

(1) Meson unparticle model 

(2) Fermion unparticle model 

(3) A free quark and anti-quark state in the Wilson fermion formalism 

The meson unparticle model and fermion unparticle model are motivated by the 
soft-wall model in AdS/CFT correspondence.^^ 

(1) The soft- wall model predicts the form of the propagator with scale dimension 
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A in the momentum space as 

iOiP)Oi-p)) = ■ (24) 

In this case, when ^ 1, ait) = 3 — 27*. This is model independent and 
universaL On the other hand, when 1, a{£) = 2 — 7* for t ^cft- This 

is model dependent. 

(2) Another plausible scenario is to treat ^p"/5'ip{x) as the non-bound state of un- 
fermions. The soft-wall model predicts the form of the propagator with scale 
dimension A/ in the momentum space as 

(vl/(p)^(-p)) = (p^7,, + m)— ■ (25) 

In this case, when mt <C 1, a{t) = 3 — 27*. On the other hand, when mt^l, 
a{t) = 1.5 — 7* for t ^ ^CFT- This is model dependent. 

Let us discuss finite sizes effects. We have performed simulations at several lattice 
sizes up to 24"^ x 96. However mn = 0.2 ^ 0.4 even for the massless quark niq ^ 
due to finite size effects. To obtain the 7* from the scaling, ea. (|5.1.2p . finite size 
effects for mn should be controlled. To estimate the finite size effects, we investigate 
them in the models. 

• We make numerical investigation of finite size effects in the case of a free 
quark and an anti-quark state in the Wilson fermion formalism (with bare 
quark mass mo = 0.20). We make a Yukawa- type fit on A^^ lattices {N = 
24, 48, 64, 128, 160, 192) using mathematica, and obtain for the exponent of the 
power and the mass of " meson" . 

p = 1.5 -f 95.8 X N-^/^ 
mjj = 0.4-0.165 X iV"i/i° . (26) 

To suppress the error of p to order 1%, we need N > 50, 000. The finite size is 
much more severe than expected. 

In the case of asymmetric lattice x 4N, there is a shoulder taking the value 
a ~ 1.5 at iVt iV, although the shape of the shoulder is not so clear as in the 
case Nf — 16. 

• In the Nf — 7 case, finite size effects are expected to be similar to those of 
mesons in the confining phase: We need N > 5 x l/mn- 

• In the Nf = 16 case, severe effects are expected, since this case is similar to the 
free case. 

12. What Kind of Theories are defined ? 

First let us consider the Nf — 16 case. The t dependency of a{t) has a shoulder 
at a 1.5, except for the case K = 0.13322 on the 24'^ x 96 lattice. We also see 
a shoulder in the Nf = 2 case at T ^ 100 Tc and T ~ 10^ Tc as well as in the free 
case. In the limit T — > 00, the quark will become a free particle state ultimately. 
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Therefore we may consider that the state at very high temperatures in the Nf = 2 
case is ahuost a free state, but is not. Then from the correspondence between Nf 
and T, which our data suggest, we conclude that our data are consistent with the 
picture that the Nf = 16 theory is close to the fermion unparticle model. 

Next let us consider the Nf — 7 case. The t dependence of a{t) is very similar 
to the Nf = 2 case at T 2Tc. We know that wave functions of meson just above 
Tc can be obtained numerically and a "meson" state is observed experimentally, 
although quarks are deconfined. This implies that the state just above Tc is a almost 
bound state, but it is not. Now, from the correspondence between Nf and T, we 
conclude that our data are consistent with the picture that the Nf — 7 theory is 
close to the meson unparticle model. 

Now we estimate 7* for Nf = 7 and Nf — 16, using the formula we have derived. 

• We assume that the Nf = 7 case is well described by the meson unparticle 
model. The value of a{t) at plateau(t = 15 ~ 31) is 0.8(1) for K = 0.1452 and 
K = 0.1459. Applying the formula a{t) = 2-7*, we have 2-7* = 0.8(1). Thus 
7* = 1.2(1). This is a crude estimate of 7* based on some assumptions. The 
estimate should be much more refined. 

• We assume that the Nf — 16 case is well described by the fermion unparticle 
model. With the value of a{t) at shoulder (i = 20 ^ 24), a{t) ^ 1.5., applying 
the formula a (i) = 1.5— 7*, we have 1.5— 7* ^ 1.5. From this we obtain 7* ^ 0.0, 
which is consistent with a small value of 2-loop calculation; 7* ^ 0.025. 

13. Solutions for Long Standing Issues at High Temperatures 

There are long standing issues at high temperatures. 

(1) The free energy of quark-gluon plasma state does not reach that of the Stcfan- 
Boltzmann ideal gas hmit even at T/Tc = 100^ 

(2) Wave functions of meson just above Tc can be obtained, although quarks are 
not confined. 

(3) What is the order of the chiral phase transition in the Nf = 2 case, 1st or 
2nd ? This is theoretically and phenomenologically important issue. Pisarski 
and WilczelP^ derived that if the order parameter of chiral symmetry satisfies 
a scaling relation with exponents of the three dimensional 0(4) sigma model, 
the order is 2nd. We showecP^ that the scaling is well satisfied using the Wilson 
fermion formalism. Although it took some time, for staggered fermion, it was 
also shown that the scaling is satisfied in the Nf = 2+1 case.l^ Therefore it 
has been believed that this issue was settled down. However, recently, Sinya 
Aoki et al. have shown that under reasonable assumptions the effect of the 
LOi(l) symmetry breaking does not contribute to some physical quantities.^i^ 
This implies that the transition is 1st order. The issue has been revised. 

Solutions for the issues are given as follows. 
(1) When r ^ Tc, quarks are not free particles, but fermion unparticlcs. Therefore 
the degree of freedom will be effectively counted less than the free fermion case. 
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(2) When T ^ Tc, "mesons" are meson unparticles. Meson unparticles become 
mesons in the hmit T ^ Tc- Therefore meson unparticles are similar to meson 
particles in some aspects. 

(3) In ref. [T71 it is assumed that physical quantities are analytic in terms of 
and TO//. However, G{t) is not analytic in terms of niq and mn, as shown in 
this article. Therefore, the assumption of analyticity does not hold. Thus the 
conclusion that the chiral phase transition in the Nf = 2 case is 2nd order is 
revised. 



14. Discussion and Conclusions 

First of all, our conjectures for the conformal theories with IR cutoff have been 
verified in the cases of Nf — 7 and Nf — 16. This fact strengthens the assertion 
that the conformal window is 7 < Nf < 16. We have further verified the Conformal 
Theories with IR cutoff in the case T/Tc > 1 for Nf — 2. We have stressed that IR 
cutoff is inherent with simulations on a lattice and QCD at high temperatures. 

The '" micro-analysis" of propagators is proposed and is shown to be very useful. 
Even in the confining phase the t dependence of a{t) contains information of RG 
flow from UV to IR. If it is affected by an IR fixed point, it will appear as remnant 
in the a{t). If it is not affected, it is difficult to estimate 7* correctly from the scaling 
rule, eq. (|5.1.2|) . 

By the '"micro-analysis" of the propagator, we observe " Nf — 7" and "T ^ Tj," 
are similar to each other, and " Nf — 16" and "T ^ Tc" are similar to each other. 
Using the correspondence between " Nf" and "T" that our data suggest, we conclude 
our data are consistent with the picture that " Nf = 7" and "T ~ Tc" are close to 
the meson unparticle model, and "iV/ = 16" and "T S> Tc" are close to the fermion 
unparticle model. 

Surely a lot of things theoretically and numerically should be done. In particular, 
the continuum limit of propagators in the large Nf case should be estimated in order 
to obtain with more confident what kind of theories are defined. To do so, we have 
to investigate finite size effects and the /3 dependence of propagators. In parallel 
with these efforts, physics implications of unparticles such as the effects in the early 
universe should be much more deepened. 
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Fig. 17. The exponent a{t) on a lattice 16^ X 64.(a) Nf = 2T ^ 2Tc.(b) Nf = 7;K = 0.1446. 



. K=0.135 
2.5 - r . mq=0.028 



win_size=3 
win_size=5 



2 

f 1.5 
1 

0.5 





3 
2.5 
2 

f 1.5 
1 

0.5 





K=0.130 
mq=0.10 



win_size=3 
win_size=5 



2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 

(a) (b) 
Fig. 18. The exponent a(t) on a lattice 16^ x64. (a) Nf = 2;T ^ 100 Tc. (b) Nf = 16; K = 0.130. 
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(a) (b) 
Fig. 19. The exponent a{t) on a lattice 16^x64. (a) Afj = 2; T ~ 10^ Tc. (b) Nf = 16; K = 0.1315. 



